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Abstract. The conjecture that the orbit-counting generating function for totally symmetric 
plane partitions can be written as an explicit product-formula, has been stated independently 
by George Andrews and David Robbins around 1983. We present a proof of this long-standing 
conjecture. 



1. Proemium 

In the historical conference Combinatoire Enumerative that took place at the end of May 1985, 
in Montreal, Richard Stanley raised some intriguing problems about the enumeration of plane 
partitions (see below) , which he later expanded into a fascinating article 9J . Most of these problems 
concerned the enumeration of "symmetry classes" of plane partitions that were discussed in more 
detail in another article of Stanley [TU]. All of the conjectures in the latter article have since been 
proved (see David Bressoud's modern classic [3]), except one, which until now resisted the efforts 
of the greatest minds in enumerative combinatorics. It concerns the proof of an explicit formula 
for the q-enumeration of totally symmetric plane partitions, conjectured, ca. 1983, independently 
by George Andrews and David Robbins ([ID]. [3] conj. 7, [3] conj. 13, and already alluded to in 
PQ). In the present article we finally turn this conjecture into a theorem. 

A plane partition 7r is an array 7r = (i^i,j)i<ij, of positive integers 7Tjj with finite sum | tt | = ^ 
which is weakly decreasing in rows and columns so that 7Tj j > 7Tj +1 j and Ttij > 7Tj A plane 
partition n is identified with its 3D Ferrers diagram which is obtained by stacking tt^ j unit cubes 
on top of the location This gives a left-, back-, and bottom-justified structure in which we 

can refer to the locations (i, j, k) of the individual unit cubes. If the diagram is invariant under the 
action of the symmetric group S3 then n is called a totally symmetric plane partition (TSPP). In 
other words, tt is called totally symmetric if whenever a location k) in the diagram is occupied 
then all its up to 5 permutations {(i,k,j),(j,i,k),(j,k,i),(k,i,j),(k,j,i)} are occupied as well. 
Such a set of cubes, i.e., all cubes to which a certain cube can be moved via S3 is called an orbit; 
the set of all orbits of n forms a partition of its diagram (see Figure [I]) . 

In 1995, John Stembridge [11] proved Ian Macdonald's conjecture that the number of totally 
symmetric plane partitions with largest part at most n, i.e., those whose 3D Ferrers diagram is 
contained in the cube [0,n] 3 , is given by the elegant product-formula 

n i+j+k-l 
i + j + k-2' 

l<i<j<k<n J 

Ten years after Stembridge's completely human-generated proof, Andrews, Peter Paule and Carsten 
Schneider [2] came up with a computer-assisted proof based on an ingenious matrix decomposition, 
but since no g-analog of their decomposition was found, their proof could not be extended to a 
proof for the g-case. A third proof of Stembridge's theorem [7], completely computer-generated, 
was recently found in the context of our investigations of the g-TSPP conjecture. We have now 
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FIGURE 1. A totally symmetric plane partition with largest part 7. The yellow 
cubes form an orbit of length 6 (corresponding to all permutations of (2,4,6)), 
the red cubes an orbit of length 3 (all permutations of (1,1,7)), and the green 
cube on location (4, 4, 4) represents an orbit of length 1 . 

succeeded in completing all the required computations for an analogous proof of the g-TSPP 
conjecture, and can therefore announce: 

Theorem 1. Let %/ S3 denote the set of orbits of a plane partition ir under the action of the 
symmetric group S3. Then the orbit- counting generating function (]2 p. 200], [101 p. 106] ) is 
given by 

1 — q i+ i +k ~ 1 



E ^ = n r 

7r<ET(n) l<i<j<k<n 

where T(n) denotes the set of totally symmetric plane partitions with largest part at most n. 



qi+j+k-2 



Proof sketch. Our proof is based on a result by Soichi Okada [H] who has shown that the theorem 
is implied by a certain — conjectured — determinant evaluation. These preliminaries are stated 
accurately in Section [2] below, followed by a description of the holonomic ansatz [14] that we are 
going to pursue. This approach relies on a kind of oracle that tells us a description of a certificate 
function c n j; the odyssey how this function has been "guessed" is described in Section [3| Once 
c„j is known, the determinant evaluation reduces further to proving the three identities ([l]), 
and |3]) stated below. In Section [4] we come full circle by proving these identities. Additional 
technical details of the proof and our computations as well as the explicit certificates are provided 
electronically on our website (http://www.risc.uni-linz.ac.at/people/ckoutsch/qtspp/). □ 

Our proof is noteworthy not only for its obvious significance in enumerative combinatorics, as it 
settles a long-standing conjecture, attempted by many people. It is noteworthy also for computa- 
tional reasons, as the computations we performed went far beyond what has been thought to be 
possible with currently known algebraic algorithms, software packages, and computer hardware. 

2. The Telemachiad 

In order to prove the g-TSPP conjecture, we exploit an elegant reduction by Okada [8] to the 
problem of evaluating a certain determinant. This determinant is also listed as Conjecture 46 in 
Christian Krattenthaler's essay [B] on the art of determinant evaluation. 



A PROOF OF GEORGE ANDREWS' AND DAVID ROBBINS' <j-TSPP CONJECTURE 



3 



Let, as usual, Si 7 - be the Kronecker delta function and let, also as usual, 



(l-^Xl-g"" 1 )-- •(!-<? 



n-k+l) 



(l-g*)(l-g*-i)-..(l- ? ) 
denote the g-binomial coefficient. Define the discrete function ctj j by: 



i + j-2 
i - 1 



i 



+ (l + q l )5 t .. J - 5 id+1 



Okada's crucial insight is that Theorem [T] holds if 



det(cijj) 



n 



l<i<j<k<n 



1 — q i +i+ k ^ 1 
1 _ gi+j+k-2 



=■■ b n (n > 1). 



So for proving the q-TSPP conjecture, it is sufficient to prove this conjectured determinant eval- 
uation. For this purpose we apply a computational approach originally proposed in |14j which is 
applicable to identities of the form 

det(aij)i<i i j< n = b n (n > 1) 

where dij and b n (with b n =/= for all n > 1) are given explicitly (as it is the case here). 

The approach rests on the following induction argument on n. For n = 1, the identity is trivial. 
Suppose the identity holds for n — 1; then the linear system 



dl,n-l 



— l,n— 1 Qn—\,n 



V o 







, n — 1 



1 / \ Cn,n J 





V) 



has a unique solution {c n ,\, ■ ■ ■ ,c n , n )- The component c n j of this solution is precisely the (n,j)- 
cofactor divided by the (n, n)-cofactor of the nxn-determinant. The division is meaningful because 
the in, n)-cofactor is just the (n— 1) x (n— l)-determinant, which by induction hypothesis is equal 
to which by general assumption is nonzero. Since the n x rt-determinant can be expressed 

in terms of the matrix entries a n j and the normalized cofactors c n j via 



b n -l I 



•n,3 n ->3 ' 



3=1 



the induction step is completed by showing that this sum evaluates to b n . 

The difficulty is that this last summation involves the function c nj (of the discrete variables n 
and j) for which we do not have an explicit expression for general n and j. In order to achieve our 
goal, we guess a suitable description of a function c n ,j and then prove that it satisfies the three 
identities 



(1) 
(2) 

(3) 



E 

3=1 



3=1 



b n -l 



(n>l), 
(1 < % < n), 

(n > 1). 



Once this is done, then by the argument given before the c raj must be precisely the normalized 
(n, j)-cofactors of the n x n-determinant and the determinant evaluation follows as a consequence. 
So in a sense, the function c n j plays the role of a certificate for the determinant identity. 
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3. The Odyssey 

In our the setting the certificate function c n j will be described implicitly by a system of linear 
recurrence equations in n and j with coefficients depending polynomially on q, q J and q n . Such 
recurrence equations can be phrased as the elements of some noncommutative operator algebras 
Q(<7, q n , qj)[S n , Sj] where the symbols S x represent the shifts x i-> x+1. If a function is annihilated 
by certain operators (viz. it satisfies certain recurrence equations), then it is also annihilated by 
all the elements in the (left) ideal generated by those operators. We speak of an annihilating ideal 
and represent such ideals by (left) Grobner bases, so that for instance ideal membership can be 
decided effectively. 

The annihilating ideals we use for representing functions are such that they uniquely determine 
the function up to some finitely many initial values. Technically, this means that the ideals 
have dimension zero and that some particular polynomial coefficients appearing in the recurrence 
system must not vanish simultaneously for infinitely many points (n,j). These recurrence systems 
are similar but somewhat simpler than g-holonomic systems |13j . which satisfy some additional 
requirements that are not needed for our proof. We have checked that all the ideals arising in our 
proof are indeed of the desired form, but in the interest of clarity we suppress a more detailed 
description of these checks here. For a complete analysis including all technical details we refer 
to the supplementary material on our website. Also in the interest of clarity, we will from now on 
identify recurrence equations with their corresponding operators. 

A priori, there is no reason why the normalized cofactors c n j should admit a recursive description 
of the kind we are aiming at, but there is also no reason why they should not. It turns out that 
they do, and this is fortunate because for functions described in this way, techniques are known 
by which the required identities ([lj, ^ and ([3| can be proven algorithmically [T3J [T2J |5J 0] . In 
order to find a recursive description for c„.j, we first computed explicitly the normalized cofactors 
c n ,j € Q(g) for a few hundred specific indices n and j by directly solving the linear system quoted 
in Section [2j Using an algorithm reminiscent of polynomial interpolation, we then constructed a 
set of recurrences compatible with the values of c„,j at the (finitely many) indices we computed. 

Polynomial interpolation applied to a finite sample u\, . . . , Uk of an infinite sequence u n will always 
deliver some polynomial p of degree at most k — 1 which matches the given data. If it turns out 
that this polynomial matches some further sequence terms Uk+i,Uk+2i ■ ■ ■ , then it is tempting to 
conjecture that u n = p(n) for all n. The more specific points n are found to match, the higher is 
the evidence in favor of this conjecture, even though without further knowledge no finite amount 
of points will ever constitute a rigorous proof. 

Very much analogously, it is possible to extract recurrence equations from some finite number 
of values c n j. The equations become trustworthy if they also hold for points (n, j) which were 
not used for their construction. In this way, we have discovered a system of potential recurrence 
equations for the c n j, which, despite being respectable in size (about 30 Megabytes), appears to 
be a rather plausible candidate for a recursive description of the normalized cofactors c n j. The 
system is available for download on our website. 

4. The Nostos 

Now we switch our point of view. We discard the definition that c n .j be the normalized (n,j)- 
cofactor and redefine the discrete function c n j as the unique solution of the guessed recurrence 
system whose (finitely many) initial values agree with the normalized (n, j)-cofactor. If we succeed 
in proving that c„j defined in this way satisfies the identities ([!]), ^ and ([3]), then we are done. 
For this purpose we provide operators which belong to the annihilating ideal of c n j or related 
ideals and have certain features which imply the desired identities. So in a sense, these operators 
play the role of certificates for the identities under consideration. 

Because of their astronomical size (up to 7 Gigabytes; equivalent to more than one million printed 
pages; corresponding to about 2.5 tons of paper), these certificates are not included explicitly in 
this article but provided only electronically on our website. Also because of their size, it was not 
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possible to construct them by simply applying the standard algorithms from [131 1121 HJ [5]. A 
detailed explanation of how exactly we found the certificates is beyond the scope of this article 
and will be given in a separate publication. But this lack of explanation does not at all affect the 
soundness of our proof, because the correctness of the certificates can be checked independently 
by simply performing ideal membership tests. Our certificates are so big that even this "simple" 
calculation is not quite trivial, but a reader with a sufficient amount of patience and programming 
expertise will be able to do it. 

We proceed by explaining the properties of the certificates provided on our website and why they 
imply 0, Q and ^. 

4.1. A certificate for ([I]). To certify that c n<n = 1 for all n > 1, we provide a recurrence in the 
annihilating ideal of c n j which is of the special form 

P7(q,q 3 ,q n )c n+7J+7 H \-p 1 (q,q : > ,q n )c n+lj + 1 +Po{q,q 3 ,q n ) = 0. 

By virtue of the substitution j i— > n, it translates into a recurrence for the diagonal sequence c n . n . 
This recurrence contains the annihilating operator S n — 1 of the constant sequence 1 as a (right) 
factor, and therefore, after checking C\ \ = c 2 .2 = • • • = c 77 = 1, it can be concluded that c„ j71 = 1 
for all n. 

Similarly, we showed that c n fl = for all n > 1 and that c n j = for all j > n. This knowledge 
greatly simplifies the following proofs of the summation identities. 

4.2. Certificates for |2]). In order to prove the first summation identity, we translate ^ into 
the equivalent formulation 



q i _ i 



i + j-2 
i - 1 



c n ,i_i - (q l + l)c ?M (l<«<n), 



taking into account that c n = 0. We provide certificates for (pf). 

To this end, let c' n i j be the summand of the sum on the left-hand side. A recursive description 
for c' n i j can be computed directly from the defining equations of c n j and the fact that the rest 
of the summand is a g-hypergeometric factor. In the corresponding operator ideal /' we were able 
to find two different recurrence equations for c' n i j which are of the special form 

P4,3<4 + 4, i+3i j H Plfl c n+l,i,j +P0,l c n,i+l,j +P0,0<^i,i,j = *n,i,j+l _ tn.i.j 

where the ^4,3, ■ ■ • ,Po,o are rational functions in < Q,{q,q l ,q n ) and t n ,i,j is a q % , q 1 , g")-linear 
combination of certain shifts of c' n i • (which are determined by the Grobner basis of V as described 
in [1]). Next observe that c' ni ^ =0 for j < (because of the g-binomial coefficient) and also for 
j > n (because c n _j = for j > n). When summing the two recurrence equations for j from — 00 
to +00, the right-hand side telescopes to and the left-hand side turns into a recurrence for the 
sum in pfj ). This method is known as creative telescoping. Finally we obtain two annihilating 
operators P\ and P2 for the left-hand side of 

For the right-hand side of ([2]]) , we can again construct an ideal of recurrences from the defining 
equations of c n j. It turns out that this ideal contains Pi and P%, so that both sides of (J2J]) are 
annihilated by these two operators. Additionally, Pi and P2 have been constructed such that 
they require only finitely many initial values to produce a uniquely determined bivariate sequence. 
Therefore it suffices to check that the two sides of ([2|]) agree at some finitely many points (n,j). 

4.3. The certificate for To certify the final identity, we rewrite ([3| equivalently into 



§) (l + 0-Cn,n-l + £ 



q n +3- l (q n +3 +q n - q- 1) 



q n - 1 



n + j - 2 
n — 1 



= (g 2ra w)n 
^ (q n ;q 2 )l 



(n > 1) 



where (a; q) n := (1 — a)(l — qa) •••(! — g" _1 a) denotes the q-Pochhammer symbol. As before, we 
use creative telescoping to provide a certified operator P which annihilates the sum on the left- 
hand side. In the present case, a single operator is sufficient because the sum depends only on a 
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single variable n (there is no i there). Using the operator P and the defining equations for c n j, we 
then construct a recurrence for the entire left-hand side, which turns out to have order twelve. As 
this recurrence is a left multiple of the second order operator annihilating the right-hand side, the 
proof of Theorem[T]is completed by checking ([3[]) for n = 1,2, ... ,12. Quod Erat Demonstrandum). 

5. Epilogue 

Paul Erdos famously (and naively) believed that every short and elegant mathematical statement 
has a short and elegant "proof from the book", and if humans tried hard enough, they would 
eventually find it. Kurt Godel, on the other hand, meta-proved that there exist many short 
statements whose shortest possible proof is very long. It is very possible that the g-TSPP theorem 
does have a yet-to-be-found proof from the book, but it is just as possible that it does not, and 
while we are sure that the present proof is not the shortest-possible (there were lots of random 
choices in designing the proof), it may well be the case that the shortest-possible proof is still very 
long, and would still require heavy-duty computer calculations. 

Be that as it may, we believe that our general approach, and our way of taming the computer 
to prove something that seemed intractable with today's hardware and software are very elegant, 
and deserve to be included in the book. 
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